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We make up-to-date QCD predictions for open charm and bottom production at RHIC in nucleon- 
nucleon collisions at -\/S = 200 GeV. We also calculate the electron spectrum resulting from heavy 
flavor decays to allow direct comparison to the data. A rigorous benchmark, including the theoretical 
uncertainties, is established against which nuclear collision data can be compared to obtain evidence 
for nuclear effects. 



Over the past few years, heavy quark production at 
colliders has received considerable attention since a va- 
riety of measurements (photon-photon, photoproduction 
and pp collisions) seemed to suggest a discrepancy, in par- 
ticular for bottom production, with respect to standard 
next-to-leading order (NLO) Quantum Chromodynamics 
(QCD) predictions. New theoretical analyses and better 
measurements have, however, largely reduced this dis- 
crepancy to the point that, in most cases, it no longer 
appears significant (see Refs. [1,2] for a review and rele- 
vant references). 

It is important to continue to validate this theoreti- 
cal framework and its phenomenological inputs 1 in new 
measurements such as the recent heavy flavor data ob- 
tained at the Relativistic Heavy Ion Collider (RHIC) by 
the PHENIX [3] and STAR [4, 5] Collaborations. Data 
taken in pp and d+Au collisions at ^/Snn = 200 GeV 
and compared to theoretical benchmark calculations will 
aid in the interpretation of heavy flavor production in 
nucleus-nucleus collisions at the same energies. In these 
heavy ion collisions, which seek to establish the existence 
of the quark-gluon plasma (QGP), a number of effects 
on heavy flavor production have been predicted. Of par- 
ticular interest are effects which modify the transverse 
momentum spectra of heavy flavor hadrons and their de- 
cay products, including energy loss [6-9], transverse mo- 
mentum broadening in both cold nuclear matter [10] and 
in passage through a hadronizing QGP [11] as well as 



collective effects such as transverse flow [12, 13]. In ad- 
dition, J/ij) regeneration in a QGP from the initial open 
charm yield has been suggested [14-16]. Thus up-to-date 
benchmark calculations of both the total charm yield and 
the transverse momentum spectra are imperative. 

The RHIC data are conveniently presented as real 
observables, either as reconstructed hadronic decays of 
charm mesons or as the heavy flavor decay electron spec- 
tra, with contributions from both charm and bottom 
hadron decays. In both cases, the transverse momentum 
spectra are presented. Such concrete observables, which 
can be directly compared to predictions of the same quan- 
tities, stand in contrast to the often adopted procedure of 
experimental 'deconvolution' to a more basic level. Such 
a deconvolution to the bare heavy quark level and the 
subsequent extrapolation to full phase space, sometimes 
involving large factors, risks biasing the data since the- 
oretical prejudice enters in both the deconvolution and 
the extrapolation, and should therefore be avoided. 

The purpose of this paper is neither to review all possi- 
ble methods to evaluate the heavy quark cross section in 
heavy ion collisions, nor to perform a detailed analysis of 
the data. Instead, besides comparing the RHIC data to 
the most up-to-date QCD predictions, we establish the 
aforementioned benchmark calculation for further com- 
parisons. To this end, we thus adhere to the rigorous 
QCD framework shown to be successful in pp collisions. 
Significant deviations from this benchmark could thus 



x We emphasize that such inputs were not chosen in order to fit the data, but rather consistently extracted from other exper- 
imental measurements. 
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signal the presence of effects specific to the high density 
environment of heavy ion collisions, such as those men- 
tioned previously. 

To make comparisons at various levels (while preferring 
the final observable), in this letter we present predictions 
of the transverse momentum, p T , distributions of charm 
and bottom quarks, the charm and bottom hadron dis- 
tributions resulting from fragmentation and, finally, the 
electrons produced in semileptonic decays of the hadrons. 
At each step, we clarify the theoretical framework as well 
as the parameters and phcnomcnological inputs. Theo- 
retical uncertainties are estimated as extensively as pos- 
sible since comparisons of data with theory should not be 
performed at the 'central value' level only but should also 
include the respective uncertainties. Our final prediction 
is thus not a single curve but rather an uncertainty band 
which has a reasonably large probability of containing 
the 'true' theoretical prediction 2 . 

The theoretical prediction of the electron spectrum in- 
cludes three main components: the p T and rapidity dis- 
tributions of the heavy quark Q in pp collisions at = 
200 GeV, calculated in perturbative QCD; fragmentation 
of the heavy quarks into heavy hadrons, Hq, described 
by phenomenological input extracted from e + e~ data; 
and the decay of Hq into electrons according to spectra 
available from other measurements. This cross section is 
schematically written as 

Ed 3 a{e) _ E Q d 3 a{Q) 
dp 3 dpQ 

® D(Q -» H Q ) ® f(H Q -> e) , 

where the symbol ® denotes a generic convolution. The 
electron decay spectrum term J(Hq — > e) also implicitly 
accounts for the proper branching ratio. 

The distribution Ed 3 a(Q)/dpq is evaluated at the 
Fixed-Order plus Ncxt-to-Leading-Log (FONLL) level, 
implemented in Ref. [17]. In addition to including the 
full fixed-order NLO result [18,19], the FONLL calcu- 
lation also resums [20] large perturbative terms propor- 
tional to a™ log (p T /m) to all orders with next-to- leading 
logarithmic (NLL) accuracy (i.e. k = n, n — 1) where to 
is the heavy quark mass. The perturbative parameters 
are the heavy quark mass and the value of the strong cou- 
pling, a s . We take m c — 1.5 GeV and mi, = 4.75 GeV 
as reference values and vary the masses over the range 
1.3 < m c < 1.7 GeV for charm and 4.5 < nib < 5 
GeV for bottom to estimate the resulting mass uncer- 
tainties. The QCD scale at five flavors, A^ 5 ), is set to 



0.226 GeV, i.e. the value provided by the CTEQ6M par- 
ton densities. The perturbative calculation also depends 
on the unphysical factorization (jUp) and renormalization 
(/j.r) scales. The sensitivity of the cross section to their 
variation can be used to estimate the perturbative un- 
certainty due to the absence of higher orders. We have 
taken \ir^f = Mo = \Jv\ + m2 as a central value and 
varied the two scales independently within a 'fiducial' re- 
gion defined by \ir.f = S,r,f^o with 0.5 < £r,f < 2 
and 0.5 < £r/£,f < 2. In practice, we use the following 
seven sets: {(&,&■)} = {(1,1), (2,2), (0.5,0.5), (1,0.5), 
(2,1), (0.5,1), (1,2)}. The envelope containing the result- 
ing curves defines the uncertainty. Finally, the uncertain- 
ties stemming from mass and scale variations are added 
in quadrature. 

These 'perturbative' inputs lead to a FONLL total cc 
cross section in pp collisions of crf e ONLL = 256^™ nh at 
VS = 200 GeV. The theoretical uncertainty is evaluated 
as described above. The corresponding NLO prediction 3 
is 244+^34 fib. Thus the two calculations are equivalent 
at the total cross section level within the large perturba- 
tive uncertainties, as expected. The total cross section 
for bottom production is cr^ ONLL = l-87±g;|? fib. 

The fragmentation functions D(c — > D) and D(b — > 
B), where D and B indicate a generic admixture of 
charm and bottom hadrons, are consistently extracted 
from e + e~ data in the context of a FONLL- type calcula- 
tion, as described in Refs. [21-23]. The charm fragmen- 
tation function [23] depends on the parameter r [24] with 
r = 0.1 for m c = 1.5 GeV, r = 0.135 for m c = 1.7 GeV, 
and r = 0.06 for m c = 1.3 GeV from e + e~ fits. Bottom 
fragmentation instead depends on the parameter a in a 
functional form by Kartvelishvili et al. [25]: a = 29.1 for 
m b = 4.75 GeV, a = 34 for m b = 5 GeV, and a = 25.6 
for nib = 4.5 GeV (see Ref. [22]). It is worth noting that 
using the Peterson et al. fragmentation function [26], 
with standard parameter choices e c ~ 0.06 ± 0.03 and 
€b ~ 0.006±0.003 would not provide a valid description of 
fragmentation in FONLL [21]. Fragmentation is numeri- 
cally performed by rescaling the quark three-momentum 
at a constant angle in the laboratory frame. This choice 
is, to some extent, arbitrary. Alternatively, one might 
rescale the transverse momentum at constant rapidity. 
While all choices are equivalent at p T ^> to, they will, in 
general, lead to different results at p T ~ to, where a large 
fraction of the RHIC data lie. The ensuing uncertainty 
is, however, not larger than the perturbative ones [22] 
and will therefore not be considered in more detail. 



2 Since a rigorous statistical estimate for this probability is not possible in the frequentist sense for theoretical errors, one 
might resort to a Bayesian analysis to quantify it. While refraining from such an effort, we roughly estimate this probability 
as ~ 80-90%. 

3 Earlier papers [10] used m c = 1.2 GeV and fj,R — /if = 1\/ P% + m 2 as reference parameters. With this choice we find 
= 427 fib. 
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FIG. 1. The theoretical uncertainty bands for the charm 
quark and D meson pr distributions in pp collisions at 
VS = 200 GeV, using BR(c —>£>) = 1. Data from the STAR 
Collaboration from d+Au collisions (scaled to pp using Abm 
= 7.5) at ^/Snn = 200 GeV, final [4] and preliminary [5], are 
also shown. 
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FIG. 2. The theoretical uncertainty bands for the bot- 
tom quark and B meson pr distributions in pp collisions at 
VS = 200 GeV, using BR(6 —>£?) = !. 



The decay of the D and B mesons into electrons is 
controlled by the experimentally measured decay spectra 
and branching ratios (BR). The spectrum for primary 
B — ► e decays has been measured recently by BaBar [28] 
and CLEO [29] . We have used a model that fits the data 
well and assume it to be valid for all bottom hadrons. 



Preliminary CLEO data on the inclusive electron spec- 
trum in semi-leptonic D decays have been shown [30]. 
We fitted this spectrum and we assume it to be identi- 
cal for all charm hadrons. Finally, the contribution of 
electrons from the secondary B decays B —> D —> e has 
also been accounted for. The relevant electron spectrum 
has been obtained as a convolution of the D — > e spec- 
trum mentioned above with a parton-model prediction 
for the b — > c decay. The resulting electron spectrum 
is very soft, suggesting that its contribution to the total 
will most likely be negligible. 

To normalize the decay spectra, we use the branch- 
ing ratios for bottom and charm hadron mixtures ap- 
propriate to this high energy regime [27]: BR(£? — > 
e) = 10.86 ± 0.35%, BR(D -> e) = 10.3 ± 1.2%, and 
BR(B -> D -> e) = 9.6 ± 0.6%. 

We first present the transverse momentum distribu- 
tions for charm quarks and charm hadrons. Figure 1 
shows the theoretical uncertainty bands for the two dis- 
tributions, obtained by summing the mass and scale un- 
certainties in quadrature 5 . Note that the band is en- 
larged at low p T due to the large value of a s at low scales 
and the increased sensitivity of the cross section to the 
charm quark mass. It is also worth noting that, due to 
the fairly hard fragmentation function, the D meson and 
c quark distributions begin to differ outside the uncer- 
tainty bands only for p T > 9 GeV. The same comparison 
is shown in Fig. 2 for bottom quarks and the subsequent 
B mesons. As a result of the harder b — ► B fragmentation 
function, the two bands partially overlap for p T ~ 20 GeV 
and beyond. 

We next consider the transverse momentum distri- 
butions of electrons from D and B decays. Figure 3 
shows the contributions from D — ► e, B — ► e and 
B — > D — > e decays as well as the total. As antici- 
pated, the B — > D — > e secondary electron spectrum is 
extremely soft, only exceeding the primary B — > e decays 
at p T < 1 GeV. It is always negligible with respect to the 
total yield. We further note that the electron spectrum 
from B decays becomes larger than that of electrons from 
D decays at p T ~ 4 GeV. The qualitative features of this 
plot are in good agreement with results obtained by the 
RHIC Collaborations using the PYTHIA event genera- 
tor [31]. 



4 Note that the experimental uncertainties on the BRs have not been included in the overall uncertainty of the predictions, 
since they are much smaller than those of perturbative origin. 

5 For example, the value of the upper curve is U(p T ) = C(p T ) + y/(M s (pr) — C(p T )) 2 + (M m (p T ) — C(p T )) 2 , where C stands 
for the central values, M s is the maximum cross section obtained by choosing m c = 1.5 GeV with the scale factors in our seven 
fiducial sets, and M m is the maximum cross section obtained with £ B = £ F = 1 and m c = 1.3, 1.5 and 1.7 GeV. 
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FIG. 3. The various components of the electron transverse 
momentum spectrum. These curves are calculated with the 
central masses and scales, i.e. m c = 1.5 GeV, m;, = 4.75 GeV 
and £,r ; f = 1. 



standard QCD calculations to heavy quark production at 
RHIC. Alternatively, a significant disagreement will sug- 
gest the need to complement this evaluation with further 
ingredients. 
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FIG. 4. The final prediction for the theoretical uncertainty 
band of the electron spectrum from charm and bottom in pp 
collisions. Preliminary data from the PHENIX [3] and STAR 
(final [4] and preliminary [5]) Collaborations are also shown. 

Finally, in Fig. 4 we show the theoretical uncer- 
tainty band for electrons coming from charm and bottom 
hadron decays at y/S = 200 GeV at RHIC. The sum of 
the three components shown in Fig. 3 corresponds to the 
central value of the band in Fig. 4. The upper and lower 
limit of the band are obtained by summing the upper and 
lower limit for each component. 

In conclusion, we have evaluated the higher order QCD 
charm and bottom quark production cross sections in 
VS = 200 GeV pp collisions at RHIC. The results are 
presented in the form of a theoretical uncertainty band 
for the transverse momentum distribution of either bare 
charm (bottom), D (B) mesons, or electrons originating 
from the decay of charm and bottom hadrons. This result 
should not be multiplied by any K factor before compar- 
ison with data. Rather, agreement within the uncertain- 
ties of the measurements will support the applicability of 



[1 
[2 
[3 
[4 
[5 
[6 

[<: 

[8 
[9 
[10 

[11 

[12 

[13 
[14 
[15 

[16 

[17 



M. Cacciari, "Rise and fall of the bottom quark produc- 
tion excess," arXiv:hep-ph/0407187. 
M. L. Mangano, "The saga of bottom production in pro- 
ton antiproton collisions," arXiv:hep-ph/0411020. 
S. Kelly [PHENIX Collaboration], J. Phys. G 30 (2004) 
S1189 [arXiv:nucl-ex/0403057]. 

J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
94 (2005) 062301 [arXiv:nucl-ex/0407006]. 

A. Tai [STAR Collaboration], J. Phys. G 30 (2004) S809 
[arXiv : nucl-ex /0404029] . 

M. Djordjevic, M. Gyulassy and S. Wicks, arXiv:hep- 
ph/0410372. 

N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. 
Rev. D 69 (2004) 114003 [arXiv:hep-ph/0312106]. 
Y. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B 519 
(2001) 199 [arXiv:hep-ph/0106202]. 

Z. w. Lin, R. Vogt and X. N. Wang, Phys. Rev. C 57 
(1998) 899 [arXiv:nucl-th/9705006]. 
R. Vogt [Hard Probe Collaboration], Int. J. Mod. Phys. 
E 12 (2003) 211 [arXiv:hep-ph/0111271]. 

B. Svetitsky and A. Uziel, Phys. Rev. D 55 (1997) 2616 
[arXiv:hep-ph/9606284]. 

V. Greco, C. M. Ko and R. Rapp, Phys. Lett. B 595 
(2004) 202 [arXiv:nucl-th/0312100]. 

Z. w. Lin and D. Molnar, Phys. Rev. C 68 (2003) 044901 
[arXiv : nucl-t h/0304045] . 

R. L. Thews, M. Schroedter and J. Rafelski, Phys. Rev. 

C 63 (2001) 054905 [arXiv:hep-ph/0007323]. 

A. Andronic, P. Braun-Munzinger, K. Redlich and 

J. Stachel, Phys. Lett. B 571 (2003) 36 [arXivmucl- 

th/0303036]. 

A. P. Kostyuk, M. I. Gorenstein, H. Stocker and 
W. Greiner, Phys. Rev. C 68 (2003) 041902 [arXiv:hep- 
ph/0305277]. 

M. Cacciari, M. Greco and P. Nason, JHEP 9805 (1998) 
007 [arXiv:hep-ph/9803400]; 

M. Cacciari, S. Frixione and P. Nason, JHEP 0103 
(2001) 006 [arXiv:hep-ph/0102134]. 



4 



[18] P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 303 

(1988) 607; 

P. Nason, S. Dawson and R. K. Ellis, Nucl. Phys. B 327 

(1989) 49 [Erratum-ibid. B 335 (1990) 260]. 

[19] W. Beenakker, W. L. van Neerven, R. Meng, 
G. A. Schuler and J. Smith, Nucl. Phys. B 351 (1991) 
507. 

[20] M. Cacciari and M. Greco, Nucl. Phys. B 421 (1994) 530 

[arXiv:hep-ph/9311260]. 
[21] M. Cacciari and P. Nason, Phys. Rev. Lett. 89 (2002) 

122003 [arXiv:hep-ph/0204025]. 
[22] M. Cacciari, S. Frixione, M. L. Mangano, P. Nason 

and G. Ridolfi, JHEP 0407 (2004) 033 [arXiv:hep- 

ph/0312132]. 

[23] M. Cacciari and P. Nason, JHEP 0309 (2003) 006 

[arXiv:hep-ph/0306212]. 
[24] E. Braaten, K. Cheung, S. Fleming and T. C. Yuan, 

Phys. Rev. D 51, 4819 (1995) [arXiv:hep-ph/9409316]. 



[25] V. G. Kartvelishvili, A. K. Likhoded and V. A. Petrov, 

Phys. Lett. B 78 (1978) 615. 
[26] C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas, 

Phys. Rev. D 27 (1983) 105. 
[27] S. Eidelman et al. [Particle Data Group Collaboration], 

Phys. Lett. B 592 (2004) 1. 
[28] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 

69 (2004) 111104 [arXiv:hep-ex/0403030]. 
[29] A. H. Mahmood et al. [CLEO Collaboration], Phys. Rev. 

D 70 (2004) 032003 [arXiv:hep-ex/0403053]. 
[30] CLEO TALK 04-42, "First Results from 

CLEO-c", J. Yelton [CLEO Collaboration], 

http://www.lns.cornell.edu/public/TALK/2004/TALK04- 

42/ 

[31] T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, 
S. Mrenna and E. Norrbin, Comput. Phys. Commun. 135 
(2001) 238 [arXiv:hep-ph/0010017]. 



5 



